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Abstract 



The Casimir energy for a conducting spherical shell of radius a 
is computed using a direct mode summation approach. An essential 
ingredient is the implementation of a recently proposed method based 
on Cauchy's theorem for an evaluation of the eigenfrequencies of the 
Q> I system. It is shown, however, that this earlier calculation uses an im- 

^ proper set of modes to describe the waves exterior to the sphere. Upon 

. making the necessary corrections and taking care to ensure that no 

00 i mathematically ill-defined expressions occur, the technique is shown 

to leave numerical results unaltered while avoiding a longstanding crit- 
icism raised against earlier calculations of the Casimir energy. 



I. Introduction 



In 1948 Casimir [1] first predicted that two infinite parallel plates in vac- 
5^ I uum would attract each other. This quite remarkable result depends only on 

the universal constants h and c and the geometry of the plate configuration, 
but is independent of such things as the electric charge. It has its origin in 
the expression for the energy E of the uncoupled electromagnetic field which 
is well known to have the form 

k,A 

where = |k| and nk,A denotes the photon occupation number in the 
mode with wave number k and polarization A. The sum is to be taken over 
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all allowed k and A. In the absence of radiation this reduces to the zero point 
vacuum energy, or 

^ = E^M'^- (1) 

k,A ^ 

This divergent expression can be rendered finite by use of a cutoff or conver- 
gence factor, and (somewhat surprisingly) made to yield measurable effects 
when boundary surfaces are used to modify the allowed set of modes which 
are to be included in (1). Thus Casimir found that there is a net attrac- 
tive force between two parallel conducting plates. This result led, naturally 
enough, to his subsequent suggestion [2] that such forces could play a signif- 
icant role in lending stability to the electron. However, the Casimir energy 
for the sphere was found to be positive, thereby implying a repulsive force 
rather than the anticipated attractive one. This result was first achieved in 
a remarkable, albeit intricate, calculation by Boyer [3] and was later verified 
by a number of authors [4-6]. 

More recently a direct mode summation approach to the problem of cal- 
culating the Casimir energy of a conducting spherical shell has been advanced 
by Nesterenko and Pirozhenko [7]. Their technique uses Cauchy's theorem 
to convert the sum over eigenmodes into an integral and was found to yield 
the same analytic expression for the Casimir energy as that obtained in refs. 
5 and 6. Since the original approach to the parallel plate geometry was also 
based on a direct mode summation, such a result provides a welcome addi- 
tion to the hterature of this subfield. On the other hand the approach of ref. 
7 has some shortcomings which prevent it from being immediately accepted 
as the natural extension of the Casimir method to the sphere. In particular 
it is shown in the following section that the standing waves exterior to the 
sphere have not been found correctly in ref. 7. After finding appropriate 
expressions for those modes the calculation is put on a more rigorous footing 
by including a cutoff function to dampen the high frequency contributions (in 
direct analogy to the Casimir calculation [1]) together with a modification of 
the contour in order to ensure the convergence of the Cauchy integral expres- 
sion. It is shown that this eliminates a spurious divergence which Candelas 
[8] asserts must be added to the result of ref. 6. The approach developed 
here applies equally well to the Dirichlet and Neumann spheres, but no ex- 
plicit calculations are carried out for these cases since earlier numerical work 
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is unaffected. An extension to (2+1) dimensional QED is considered in a 
subsequent section with the result that there is an intrinsic divergence in the 
Casimir energy in this case. 



II. Mode Summation for a Spherical Shell Using Cauchy's 

Theorem 

In evaluating the Casimir energy of a sphere it should be noted that each 
mode is 2Z + 1 fold degenerate. Thus (1) becomes 



OO -1 CX) 

1=1 ^ n=l A=l,2 



(2) 



where the eigcnfrcquencies uj^^^ are determined by imposing appropriate 
boundary conditions on the multipole fields. In particular, for the case of a 
spherical shell of radius a, the transverse electric (i.e., A = 1) modes are 



ji{uja) = 



(3) 



Aiji{ua) + Bini{uja) = 
and the transverse magnetic (i.e., A = 2) modes are 

d 



dr 







dr 



r[Ciji{ujr) + Dimiujr)] 







(4) 



(5) 



(6) 



where ji and rii are the spherical Bessel functions. The constants Ai, Bi, 
Ci, and Di (more specifically, the ratios Bi/Ai, Di/Ci) are determined by 
prescribing the correct asymptotic behavior at large r. It should be noted 
that Eqs. (3) and (5) determine the interior (r < a) modes while (4) and (6) 
specify the exterior (r > a) modes. 

The coefficients which appear in the exterior mode equations are deter- 
mined by enclosing the entire system within a second concentric conducting 
sphere of radius R. It is straightforward to verify that Eqs. (4) and (6) can 
then be written as 
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and 



ji{ijja) + tan6ini{u!a) = 



— r[ji{ujr) - cot5ini{ur)] 
ar L 

respectively, where for sufficiently large R 



r=a 



= 




(7) 



This is to be contrasted with the exterior mode solutions of ref. 7 where Eqs. 
(4) and (6) have been replaced by those of Stratton [9], namely 



While such boundary conditions are appropriate to determine the complex 
frequencies for a radiating sphere, they are not valid to specify the necessarily 
real frequencies which contribute to the vacuum energy [10]. 

Having formulated the conditions for the eigenfrequencies of the system 
it now remains to be shown how Cauchy's theorem can be applied to the 
evaluation of the Casimir energy. In analogy to ref. 7, the eigenfrequency 
equations are defined as 



h\''{uja)^0 



and 




where h, ' denotes the Hankel function 



hi '{z) = ji{z) + ini{z). 




ji{z) +tsinSi{z)ni{z) 



z[ji{z) - cot5i{z)ni{z)] 
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where z = ua and 

Si{z) = z{R/a) - ^. 
Another useful definition is what might be termed the f-product, or 

a 

where a factor of z"^ has been included for convenience in order to make fi{z) 
finite at z = 0. Defining 2; to be a complex variable, fi{z) is seen to be an 
analytic function of z, the zeros of which correspond to the eigenfrequencies 
of the system [11]. 

It follows from Cauchy's theorem that for two functions fi{z) and ^(z) an- 
alytic within a closed contour C in which fi{z) has isolated zeros at xi,X2, ■■■Xn, 
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Choosing (/){z) — ze~"'^ where cr is a real positive constant thus leads to 

Upon combining Eq.(8) with Eq.(2) the Casimir energy becomes 

1 °° 1 /■ 
= lim ya + -) i dze-^'z—lnfiiz). (9) 

It should be noted that the factor of e~"^ plays the role of a cutoff function 
which effectively suppresses the high frequency contributions to the Casimir 
energy. While such a cutoff function is generally invoked in calculations based 
on parallel plate geometry, it was not employed in refs. 6 and 7. It will be 
seen, however, that it plays an important role in ensuring that all integrals 
encountered in Eq. (9) are well defined. 

Before the calculation of the Casimir energy can be completed, it is nec- 
essary to specify a contour for the integration. An appropriate contour C 
(displayed in Fig.l) can be conveniently broken into three parts. These con- 
sist of a circular segment Ca and two straight line segments Fi and Since 
the r contours are oriented at a nonzero angle with respect to the imag- 
inary axis, it follows that the contribution to Ca is, by virtue of the cutoff 
factor, bounded by exp{—aAsin4>) where A is the radius of the circular arc. 
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Figure 1: Contour integral in complex plane. 



Since the logarithm in the Casimir energy expression grows at most alge- 
braically, it follows that the contribution to Ca vanishes exponentially in the 
limit of large A provided that 7^ 0. 

Along Fi and F2 the quantity tanSi becomes i and —i respectively for 
sufficiently large R. Thus on Fi, where y — \z\, 

with the corresponding terms for F2 obtained by complex conjugation. This 
allows fi{z) to be written, up to an overall normalization, in terms of the 
modified Bessel functions as 



= I,{ye-"f>)K,{ye-"^) 

ri 
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where v — I -\- \. Upon noting that the contributions along Fi and r2 are 
complex conjugates of each other, the Casimir energy becomes 

Ec^ + 0)^66-^^ / dyexp{-iaye-"^)y-lnfi{ye-''^). 

na j~y I Jo ay 

Straightforward manipulation [5] of the argument of the logarithm and rescal- 
ing of the integration variable allows this to be recast in the form 



Ec = -— + a^'^^'^'Vo dyexp{-ivaye-'t')y-ln^ - [\i{yve-"f)f 

(10) 

where 

dy 

Upon using the uniform expansions of lyii'y) and Ky^vy) for large u [12] it 
is found that the integral over y is finite for cr = 0. Adding and subtracting 
the asymptotic form of the integrand for large the Casimir energy assumes 
the form 

Ec = Efin + E^, 

where the finite part is 



Efin 



00 



E I dy [(/ + -)Hn[l - {\{uy)f\ + -(1 + y^) 
and the cutoff dependent part is 



2\-3 



1 00 /-OO 

= -^11^^^""^ dyexp{-waye-''^)y—{l + y^e''^''^)-^ 
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In arriving at the independent form for Efm one has in succession taken the 
hmit cr — > and replaced y with ye"^"^. The resulting integral over the interval 
< y < e~*'^oo is then trivially seen to be equivalent to the corresponding 
integral over the real interval < y < oo. The quantity Ef^n has been 
evaluated numerically in ref. [6] and consequently no attempt is made here 
to obtain that result independently. 

Although the integral which appears in E„ is finite even for a = 0, the 
limit cannot be taken at this point since it would yield the divergent result 

o oo -1 

Consequently, E^^ is to be evaluated for finite a before taking the limit of 
vanishing cutoff. Since this is clearly a delicate limiting process, some care 
is warranted. It should first of all be noted that the cutoff factor used here 
was initially introduced as a real exponential prior to its emergence in E„ 
as the term exp{—ii'aye~^^) . Thus there is a damping term of the form 
exp{—i'aysm0). This may be contrasted with ref. [6] which objects to using 
an explicit cutoff but finds, by using a temporal point separation, that an 
oscillatory factor of the form e*'^''^ can appear in E„. It is significant that 
such a term (unlike that found here) has no residual exponential damping 
effect. This has led to a criticism of the ref. 6 result by Candelas [8] who 
asserts that because the waves of different / all contribute an amount of equal 
absolute value to the Casimir energy, an explicit cutoff must be introduced 
in the sum over /. His result is that such an insertion adds a divergent term 
to the Casimir energy that is proportional to /q°° dv (where is a frequency). 
Since, however, the calculation presented here has included a cutoff ah initio 
there should be no need for the inclusion of an additional cutoff of the type 
used in [8] . This can now be verified by explicitly carrying out the summation 
over / to give 



27ra Jo (1 + y2g-2i<^j4 exp{iaye-'-'t>) - 1 ' 

For nonzero the sum is well defined and thus requires no additional cutoff 
factor. Upon expanding the cr dependent terms and dropping those which 
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go to zero as some positive power of a the above expression reduces to 

27ra Jo ^(l + y2e-2j<A)4L (7ye-*<^^ 

As before, let y e^'^y and perform a rotation of the resulting contour. This 
allows the term singular in cr to be dropped (it is purely imaginary) and 
leaves the finite result ^ 

„=o 64a ' 

Upon combining this with the numerical result of ref. [6] for Efin, the ex- 
pression for the Casimir energy is obtained, namely 

Ec = 0.09235/2a. 

As a final comment concerning the conducting sphere, it should be noted 
that the method presented here can equally well be applied to the Dirichlet 
and Neumann spheres. The results for these cases can be rigorously obtained 
using the technique of mode summation in conjunction with Cauchy's theo- 
rem. As it would consist merely of a rewriting of ref. [7] using the approach 
developed here, this is left as an exercise for the interested reader. 

III. Casimir Effect of the Electromagnetic Circle 

A natural question which arises subsequent to a consideration of the elec- 
tromagnetic Casimir effect for a sphere is whether similar results can be ob- 
tained in higher and lower dimensions. As will be shown here this question 
can readily be answered at least in the latter case. Although the electromag- 
netic field is trivial (i.e., there is no photon) in one spatial dimension, for the 
two dimensional case a procedure analogous to that of the preceding section 
can be carried out. This requires the extraction of appropriate boundary 
conditions from the field equations. The latter are contained in the covariant 
equations 

a^F^^ = (12) 
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where J'^ = {Ji, p) with i = 1,2. From (11) it follows that 

S = V X A 

and 

E = -doA - VA° 

where B is the (scalar) magnetic field and E is the (two- vector) electric field. 
Correspondingly, Eq.(12) yields 

V-E = p 

and 

-^oE + WB^3 

where (V)j = ejjVj. 

In the case that all fields have time dependence e"*"^* these equations 
imply that in current free regions 

{V^ + 00^)3 = 0. 

Thus the B field can be written in the form 

oo 
m=—QO 

with the corresponding E field given by 

-iujF. = VB. 

Since 

V X E = -doB, 

it follows that for a conducting circle of radius a the appropriate boundary 
condition is r x E|r=o — 0, or equivalently, ^B\r=a — 0. In other words, the 
Casimir effect in (2+1) dimensional QED is equivalent to that of a scalar 
field satisfying Neumann boundary conditions. 

For the inside modes of the system it follows that the eigenfrequencies 
are determined by the condition 



= 0. 

r=a 



10 



The outside modes require that in analogy to the three dimensional case the 
system be enclosed in a large circle of radius R which is subsequently taken 
to infinity. Thus the boundary condition in the exterior region takes the form 







where Si is given by Eq.(7). Upon repeating the procedure of the preceding 
section it is found that the Casimir energy of the circle has the form 

I °° r°° d 

Ec = -7r- E Ree-M dyexp{-iaye-^^)y-ln[yI'^{ye''^)K'^{ye-'% 
2na Jo dy 



m=—oo 



The question which needs to be addressed is whether upon performing the 
sum over m a finite result for the Casimir energy can be obtained in the 
cr — limit. Using the uniform expansions for the modified Bessel functions 
for large \m\ and considering only those m values with absolute value greater 
than some large integer M, it is found that 



2 oo 

Ec ^ — me' 

m=M 



-ma 



This yields the divergent result 

in the limit of vanishing a, a result at variance with the claim in ref. [13] 
of a finite Casimir energy. On the other hand, it is not unexpected in view 
of the divergence found by Bender and Milton [14] for the two dimensional 
Dirichlet Casimir effect. 
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IV. Conclusion 



Since the inception of this subfield, there have been a number of ap- 
proaches developed for calculating the Casimir energy of a conducting sphere 
in vacuum. The first successful evaluation consisted of analytically summing 
the individual mode contributions. This method, closely modelled after the 
much simpler parallel plate problem, was the most obvious at the time but 
included tedious algebraic manipulations. More advanced methods followed 
which rendered the involved mode summation approach obsolete prior to its 
re-examination in ref. 7. The introduction of Cauchy's theorem (replacing 
the mode summation with an integration in the complex frequency plane) 
opened the door to a revitalization of this method. However, the formalism 
was based upon an improper treatment of the external modes which led to 
their being complex and finite in number rather than real and infinite as re- 
quired. A correct set of equations for the exterior modes was obtained upon 
placing the entire system in a second conducting sphere whose radius R was 
eventually allowed to go to infinity. 

A second problem associated with the Cauchy integral method consisted 
in the fact that the neglect of the semicircular part of the contour integral was 
problematical. This difficulty was eliminated by invoking a cutoff function 
(as in the original parallel plate calculation), and by introducing simultane- 
ously a kink into the integration contour. This allowed the integral to be 
evaluated along the imaginary axis. It also led to a mathematically well de- 
fined summation over partial waves which eliminated the need for a cutoff of 
the type proposed by Candelas. Thus an unobservable but divergent contri- 
bution to the Casimir energy which he claimed to be essential was found not 
to be relevant. It was noted that the method developed here applies equally 
well to the Dirichlet and Neumann spheres, although no explicit calculations 
have been carried out. Finally, an extension to the case of (2+1) dimen- 
sional QED was noted with divergent results being obtained in the limit of 
vanishing cutoff. 
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